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ABSTRACT 

The Cepheus B (Cep B) molecular cloud and a portion of the nearby 
Cep 0B3b OB association, one of the most active regions of star formation 
within 1 kpc, have been observed with the IRAC detector on board the Spitzer 
Space Telescope. The goals are to study protoplanetary disk evolution and pro- 
cesses of sequential triggered star formation in the region. Out of ~ 400 pre-main 
sequence (PMS) stars selected with an earlier Chandra X-ray Observatory obser- 
vation, ~ 95% are identified with mid-infrared sources and most of these are 
classified as diskless or disk-bearing stars. The discovery of the additional > 200 
IR-excess low- mass members gives a combined Chandra+ Spitzer PMS sample 
that is almost complete down to 0.5 Mq outside of the cloud, and somewhat 
above 1 Mq in the cloud. 

Analyses of the nearly disk- unbiased combined Chandra and Spitzer selected 
stellar sample give several results. Our major finding is a spatio-temporal gra- 
dient of young stars from the hot molecular core towards the primary ionizing 
O star HD 217086. This strongly supports the radiation driven implosion (RDI) 
model of triggered star formation in the region. The empirical estimate for the 
shock velocity of ~ 1 km/s is very similar to theoretical models of RDI in shocked 



^Department of Astronomy & Astrophysics, 525 Davey Laboratory, Pennsylvania State University, Uni- 
versity Park PA 16802 

^Center for Exoplanets and Habitable Worlds, 525 Davey Laboratory, Pennsylvania State University, 
University Park PA 16802 

•^Max-Planck-Institut fur Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany 

''Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 



- 2 - 



molecular clouds. The initial mass function (IMF) of the lightly obscured trig- 
gered population exhibits a standard Galactic field IMF shape. The unusually 
high apparent value of > 70% star formation efficiency inferred from the ratio 
of star mass to current molecular gas mass indicates that most of the Cep B 
molecular cloud has been already ablated or transformed to stars. Contrary to 
the current RDI simulations, our findings indicate that star formation triggering 
by HII region shocks is not restricted to a single episode but can continue for 
millions of years. Other results include: 1. agreement of the disk fractions, their 
mass dependency, and fractions of transition disks with other clusters; 2. con- 
firmation of the youthfulness of the embedded Cep B cluster; 3. confirmation of 
the effect of suppression of time-integrated X-ray emission in disk-bearing versus 
diskless systems. 

Subject headings: ISM: individual (Cepheus B cloud) - open clusters and associ- 
ations: individual (Cepheus 0B3) - stars: protoplanetary disks - stars: pre-main 
sequence - stars: formation - X-rays: stars 



INTRODUCTION 



Cepheus B (Cep B) is a molecular core located at the edge of the Cepheus gian t molec- 



ular cloud at a dista nce around 725 pc and lying 2.6° above the Galactic Plane fjSargent 



19771 : IYu et al.lll996l ). A handful of en ibedded young s t ars were found in Cep B from radio 
continuum and infrared (IR) studies (IFelli et al.lll978l : iTesti et al.lll995l ). The unobscured 
stellar OB association Cep 0B3 lies arou nd Cep B; the younger subgrou p, Cep 0B3b, lies 
closest to the cloud as shown in Figured] (jBlaauwlll964l : iKun et al.ll2008l ). For many years, 
the Cep 0B3 association has been considere d to be a good exaraple of large-scale sequential 
star formation in accord with the model of lElmegreen fc Ladal (119771 ) where stellar winds 
and supernova remnants of an older stellar cluster compress and trigger a second generation 
of star formation in nearby molecular cloud cores (jSargentlll979l ). 



The interface between the molecular cloud and the Cep 0B3b star association is clearly 
delineated by the optically bright HII region Sharpless 155 (S 155), where cloud material is 
ionized and heated by the radia tion field of the 07 star HP 217086, Bl star HD 217061 and 
perhaps other cluster members (jPanagia fc Thunj|l98ll : lBeuther et al.ll2000l ). Figured] shows 
the spatial relationship of the cloud, HII region, and exciting stars. Unlike in the Orion 
Nebula where the HII region lies between the cloud and ourselves, the photodissociation 
region at S 155 is favorably oriented to reveal the progression of star formation. Following 
the triggered star formation model, we expect the surface of the cloud to be eroded by the 
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early-type stars so that the cloud edge moves eastward across the observer's field of view 
with new stars emerging from the obscuring molecular cloud. Sources located within but 
near the edge of the molecular cloud would represent a new generation of star formation 
triggered by the HII region shock propagating into the cloud. 

This scenario of triggered star formation has been recently strengthened by the discovery 
using the Chandra X-ray Observatory of > 300 lightly-obscured low-mass pre-main sequence 
(PMS) members of the Cep 0B3b cluster located outsid e of the cloud, and a rich population 
of ~ 60 PMS stars in the cluster embedded within Cep B (iGetman et al.ll2006l ). X-ray surveys 
are particularly effective in discriminating disk-free PMS populations from Galactic field stars 
wh ich often badly conta minate infrared (IR) surveys of young stellar clust ers (see review 



bv iFeigelson et al.l 120071 ). Using 2MASS counterparts of the X-ray sources, iGetman et al. 
( 120061 ) found that PMS stars in the embedded cluster are more likely (26% vs. 4%) to 



have K-hand excesses from heated inner protoplanetary disks than stars in the unobscured 
Cep 0B3 region. This supports both, youthfulness of the embedded Cep B population and 
the prevalence of planet forming disks in the embedded cluster. 

We seek here to elucidate the relationships of disks and environments in this region using 

mid-IR photometry from the Spizer Space Telescope which is exceptionally well-suited for de- 

tecting disks around low-m ass members of young clusters (e.g. Allen et al. 112 004 : Sicilia-Aguilar et al 
20061 : iLuhman et al.ll2008bl . and references therein). Measures of disk fractions as functions of 
stellar mass and age in different star-forming environments can provide better understanding 
of the evolution of circumstellar disks in different environments and thus have direct astro- 
physical importance in evaluating the conditions for planet formation. X-ray PMS selection is 
very complementary to IR surveys. As X-ray emission from PMS stars is based on enhanced 
solar-type magnetic reconnection events rather than disk or accretion proces ses. X-ray selec- 



tion delivers rich and clean samples o f diskless stars missed by IR selection (IFeigelson et al. 



20071 ). The X-ray stars identified by IGetman et al.l (120061 ) in combination with discovered 
here Spitzer disk-bearing stars constitute a valuable sample to measure protoplanetary disk 
properties in different radiative environments and to study the triggering process inside the 
Cep B molecular cloud. We use the spatial distribution of disks around young stars of the 
selected Chandra and Spitzer sample to test concepts of cloud ablation and the radiative 
driven implos ion (RDI) model of triggered star formation on the edges of HII regions (e.g., 
Bertoldilll989l ). 



In ^ we report the Spitzer observations using the InfraRed Array Camera (IRAC) 
detector. Classification on disk-bearing and diskless X-ray emitting PMS stars is provided 
in ^and the discovery of the new population of non-Chandra IR-excess members is given 
in ^ The IMF and mass completeness limits of the combined Chandra and Spizer stellar 



-4- 



samples of different parts of the region are considered in ^ Tlie results on the disk evolution 
of the selected PMS sample are presented in ^ We end in ^ with the implications of the 
new observational findings for the star formation process in the region. We adopt a distance 
of 725 pc to the star forming complex and infer an age of 2 — 3 Myr for the Cep 0B3b stars 
and an age between < 1 and 2 Myr for the embedded Cep B stars, as discussed in Appendix 
Rl Appendix [B] presents the relationship between X-ray activity, stellar mass and disks; 
these results confirm those found in other young stellar clusters. 



2. CHANDRA AND SPITZER OBSERVATIONS 



2.1. X-ray Data 



T he X-ray observations of the Cep B / 0B3b and their data analysis are described in de- 
tail by lGetman et al.l (120061 ). The 30 ks exposure was obtained on 2003 March 11.51 — 1 1.88 
with the Advanced CCD Imaging Spectrom eter (ACTS) detector (IGarmire et al.ll2003l ) on- 
board NASA's Chandra X-ray Observatory (IWeisskopf et al.ll2002l ) as part of the ACIS In- 
strument Team's Guaranteed Time Observions (Obsid #3502, G. Garmire, PI). The imaging 
array (ACIS-I) consists of four abutted 1024 x 1024 pixel front-side illuminated charge- 
coupled devices (CCDs) covering a bout 17^ x 17' o n the sky. Following data reduction based 
on the ACIS Extract IDL script, iGetman et al.l (120061 ) obtained an X-ray catalog of 431 
sources, most with < 0.5" positional accuracy (Figure [2]). Using 2MASS counterparts, 89% 
of the X-ray sources are confidently associated with cluster PMS members of the region. The 
remaining X-ray sources are mostly extragalactic active galactic nuclei. 



2.2. IRAC Data 

The mid-IR observat ion of Cep B and Cep 0B3b was obtained on 2007 February 18 



with the IRAC detector (iFazio et al.l |2004| ) on NASA's Spitzer Space Telescope in the 3.6, 
4.5, 5.8, and 8.0 fim channels. This is a General Observer project (program identification 
#30361; J. Wang, PI). Two adjacent fields subtending ~ 20' x 15' in channel pairs 3.6/5.8 fim 
and 4.5/8.0 /xm were centered, respectively, on a = 22^57"09!7, 6 = +62°37'52'.'3 (J2000) 
and a = 22^56'^20:3, 6 = +62°41'43'.'0 (J2000). The orientation of the fields long axis is 
P. A. = 36.3°. These fields were imaged with a 3 x 4 mosaic of adjacent positions separated 
by 290", and images were taken at 12 dithered positions for each of the 12 cells of the mosaic. 
The images were taken in the high-dynamic-range mode with long effective exposures of 96.8 s 
in 3.6, 4.5, 5.8 /xm channels, and two 46.8 s exposures in the 8.0 /im channel, as well as short 
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effective exposures of 0.4 s in all 4 channels. Full coverage in all four channels is obtained 
for a ~ 13' X 15' area, which in turn covers ~ 65% of the 17' x 17' Chandra's ACIS-I field, 
as shown in Figures [T] and [2J 

Basic Calibrated Data (BCD) products from the Spitzer Science Center's IRAC pipeline 
version S15.3.0 were automatically treated with the WCSmosaic IDL package developed by R. 
Gutermuth for the IRAC instrumental team. Starting with BCD data products, the package 
mosaics individual exposures while treating bright source artifact s, cosmic ray rejection , 



distortion correction, subpixel offsetting, and background matching (IGutermuth et al.ll2008l ). 
We selected a plate scale of 0.86" for the reduced IRAC mosaics, which is the native scale 
divided by \/2. 

Aperture photometry of IRAC sources was obtained using the IRAF task PHOT. For 
the most of the X-ray and non-Xray selected sources, an aperture radius of 4 pixels 
(3.44") was used with an adjoining sky annulus width of 1 and 6 pixels (0.86"and 5.16") for 
the 3.6/4.5 and 5.8/8.0 fim bands, respectively. For the crowded sources, aperture radii of 2 
or 3 pixels (1.72"and 2.58") and adjoining sky annulus width of 1 pixel (0.86") in all 4 bands 
were used. We adopted zero point magnitudes (ZP) of 19.670, 18.921, 16.855 , and 17.394 



in th e 3.6, 4.5, 5.8, and 8.0 fim bands, where M = —2.5log{DN/ sec) + ZP (IReach et al. 



20051 ). The total aperture corrections applied to our measurements are: 0.185, 0.175, 0.165, 
0.250 in the case of 4 pixel aperture; 0.275, 0.305, 0.515, 0.735 in the case of 3 pixel aperture; 
and 0.655, 0.735, 1.005, 1.010 in the case of 2 pixel aperture in the 3.6, 4.5, 5.8, and 8.0 fim 
bands, respectively. The reported photometric errors include Poisson errors i n the source 



and background emission plus a 2% uncertainty in the calibration of IRAC (IReach et al. 



20051). 



2.3. IRAC Counterparts to X-ray Sources 

The initial list of point sources in these images was produced by running IRAF task 
STARFIND with a low detection threshold in order to identify as many sources as possible, 
permitting false detections in this stage of the analysi s . An automated cross-correlation 



between the Chandra source positions of iGetman et al.l (120061 ) and IRAC candidate source 



positions was made using a search radius of 2" within ~ 6' of the ACIS field center, and 
a search radius of 3.5" in the outer regions of the ACIS field where X-ray source positions 
are more uncertain due to the deterioration of the Chandra telescope point spread function 
(Figure [2]). This was followed by a careful visual examination of each source in both bands 
to remove dubious sources and associations. 
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Of the 431 X-ray sources within the 17' x 17' Chandra ACIS-I field flGetman et al.ll2006r ). 
396 (92%) he within the 3.6/5.8 jjm IRAC coverage and 412 (96%) he within the 4.5/8.0 /im 
IRAC coverage. Out of [396, 412, 396, 412] X-ray sources in the [3.6, 4.5, 5.8, 8.0 fim] 
channels, [384, 399, 369, 385] are identified with IRAC sources. The IRAC detection rate 
of X-ray sources is thus ~ 97% for 3.6 and 4.5 /xm channels and ~ 93% for 5.8 and 8.0 /im 
channels. 



2MASS associations are obtained as described in lGetman et al.l (12006 ). Opt ical coun 



terparts were identified from the Vic photometric catalog of iMayne et al. 



(I200JEI. Unlike 



the search for the IR counterparts, we lacked the ability to visually examine the Chandra- 
optical matches. We therefore compared the C/ianrfr a-optical offset (Offi in Tabled]) with 
the C handra-2M AS S and Chandra-IKAC positional offsets (Off2 and Offs). Accepted op- 
tical matches typically have |Offi — Off2,3| ^ 0.2". Out of the total of 431 Chandra X-ray 
sources, 368 are in the field of view of the optical survey and 321 of them have optical 
counterparts. The detection rate of the X-ray sources in the optical survey is thus ~ 87%. 

There are [12, 13, 27, 27] X-ray sources in the [3.6, 4.5, 5.8, 8.0 /im] channels lacking 
mid-IR counterparts. Roughly half of these lie in the interface region between the Cep B 
molecular cloud and S 155 HII region where the IRAC point source sensitivity is reduced 
by the high background nebular emission from heated dust (Figured]). These are probably 
true cluster members, but are omitted from the science analysis below. The other half of 
these unidentified X-ray sources is uniformly distributed outside of the Cep B cloud with 
weak X-ray fiuxes and high X-ray median energ ies consistent with ex tragalactic background 
sources. Most of these are listed in Table 4 of iGetman et al.l (120061 ) as likely extragalactic 
contaminants. 

The results of this search for steUar counterparts of the Chandra sources are presented in 
Tabled! Listed by Chandra sequence number, the table gives offsets between the optical/IR 
stars and the X-ray sources, VIcJHKg and IRAC photometry, and photometric flags. The 
flnal column gives the aperture used in the photometry for each IRAC band; a '9' value 
indicates that the X-ray source lies outside the IRAC coverage in that band or the IRAC 
source is too weak or strongly contaminated with no reliable photometry. 



^The optical survey encompasses most of the Chandra- ACIS-I field and extends further from the ACIS-I 
^ 15' to the north-west and ^ 20' to the west. 
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3. PRE-MAIN SEQUENCE DISK CLASSIFICATION 



We base the evolutionary classification of X-ray emitting PMS stars in the Cepheus 
region on a comparison of their IR spectral energy distributions (SEDs) with t he SEDs of 



PMS stars in the well-studied IC 348 cluster in the Perseus molecular cloud (ILada et al. 



20061 ) ■ IC 348 stars have roughly the same 2 — 3 Myr age as Cep 0B3b stars. The categories 
are: Class I (protostar, disk and envelope), Class II (PMS star and accretion disk), and 
Class III (PMS star with weak or absent disks). A simple characterization of IR SEDs 
involves the SED spectral slope a of the SED a = d\og{\Fx) / d\og{\) over the range 2 < 
A < 20/im. Class I sources have < a < 3, Class II sources have — 2 < a < 0, a nd Class III 
sources have — 3 < a < —2. Using optical and IR data for the IC 348 PMS stars, iLada et al. 



(120061 ) showed the utility of the SED slope measured in the IRAC wavelength range from 
3.6 to 8.0 fim. Comparing op tical and IRAC SEDs to disk models calculated from a Monte 
Carlo radiative transfer code (IWood et al.ll2002l : IWalker et al.l 120041 ). they find that Class II 
systems with optically thick disks have dereddened IRAC slopes aa > —1.8, "transition disk" 
(TD) systems with inner disk holes or optica lly thin inne r disks have —2.56 < < —1.8, and 
diskless Class III systems have ad < —2.56. iLada et al.l (120061 ) also provide useful empirical 
templates of median SEDs of IC 348 PMS stars for each evolutionary class. 

We proceed with the disk classifications through the comparison of the observed Cepheus 
SEDs in 2MASS+IRAC IR bands to the (de)reddened median SED templates of IC 348 PMS 
stellar photospheres. We further confirm these classifications through calculations of the SED 
slope a, and locations in IRAC color-color diagrams. Results of our disk classification are 
shown in Table [2j 



3.1. Comparison with IC 348 



First, a rough estimate of individual stellar masses of the Cepheus PMS stars is made 
by de-reddening star positions in the J versus J — H color-magnitude diagram to the 3 Myr 
and 2 Myr theoretical is ochrons for the Ce p B3b/S 155 and Cep B subregions, respectively, 
based on the models of ISiess et al.l (120001 ) and iBaraffe et al.l (Il998l . see Appendix A). These 
individual stellar photometric mass and extinction estimates are listed in Table [H Although 
individual mass estimates may be subject to significant uncertainties and may not be always 
compatible with masses obtained by more accurate methods such as optical s pectroscopy, 



this approach should be adequate to assign an IC 348 SED template (Table 4 of iLada et al. 
20061 ) to individual Cepheus PMS stars. 



Second, it is important to note that the IC 348 template SEDs are observed (not de- 
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reddened) and exhibit an extinction of Ay ~ 2 — 2.5 mag (iLada et al.ll2006l ). Thus for each 
of the Cepheus stars, an IC 348 template SED at the corresponding spectral sub-class is first 
reddened and de-reddened within the wide but finit e range of extinction s using reddening 
relationships described in Appendix A for JHKg and lFlaherty et al.l (120071 ) for IRAC bands. 
The (de)reddened IC 348 template which, after normalization to an observed J-band Cepheus 
SED point, is the closest to an observed if-band Cepheus point (both J and if-band points 
are considered to be indicators of a pure photospheric emission) is recognized here as the 
best-fit to Cepheus data. Third, the best-fit IC 348 (de)reddened median template is visually 
matched to the observed Cepheus SED, allowing classification of the Cepheus stars. 

Out of the total of 431 X-ray Cepheus sources, our classification procedure yields 215 
diskless Class III systems (labelled "NoD" in Table [2| and 139 Class II or I disk-bearing sys- 
tems ("DSK"). Based on the membership analysis of iGetman et al.l (120061 ). 24 X-ray sources 
without IR counterparts are possible extragalactic contaminants ("EXG") and 13 are fore- 
ground candidates ("PRC"). Forty Chandra sources with IR detections which are unreliable 
or insufficient for classification have an uncertain classification ("UNC"). We further subjec- 
tively classify 20 stars as Class II/III transition disks (labelled as "TD" and "DSK"); these 
follow the photospheric SED at wavelengths shorter than at least ~ 3.6/im. Four Chandra 
sources ( 314, 322, 328, and 390) embedded in the Cep B core have extreme SED excesses 
with a > and are classified as Class I stars. 

Examples of diskfree, disk-bearing, and transition SEDs are shown in Figures The 
complete set of SEDs for 354 X-ray emitting Cepheus PMS stars classified here is available 
on-line as a Figure Set. The on-line atlas contains 161 diskfree stars of the Cep 0B3b region 
(labelled as "NoD" and "CepOB3b" in Table [2]); 75 disk-bearing stars of the Cep 0B3b 
("DSK" and "CepOB3b"); 36 diskfree stars of the S 155 region ("NoD" and "S 155"); 26 disk- 
bearing stars of the S 155 region ("DSK" and "S155"); 18 diskfree stars in the Cep B region 
("NoD" and "CepB"); and 38 disk-bearing stars in the Cep B region ("DSK" and "CepB"). 
Quantitative information derived from the fitting process is given in Table [2] including the 
masses and absorptions estimated from the J — {J — H) color-magnitude diagram (columns 
4 and 5) and the SED classification obtained by comparison with IC 348 templates (columns 
6 and 7). 



3.2. Calculation of SED slope a 

We performed least-squares linear fits to the log(AFA) values in the four IRAC wave- 
length band to obtain the observed (not de- reddened) SED spectral index, oq. These are 
tabulated in column 2 of Table [2J Figure [H] shows the histogram of 354 oq values for the X-ray 
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selected Cepheus PMS stars compared to 299 IC 348 PMS stars from lLada et al.l (120061 ) 



Two important features are noted. First, the Cepheus histogram shows a bimodal 
distribution with peaks at ao ~ —1.3 and —2.7 corresponding to disk-bearing and diskless 
star samples, respectively. The IC 348 histogram ha s the same peaks, and similar bi-modality 
i s seen in the PMS pop ulations of the Cha I cloud (ILuhman et al.ll2008al ) and a Ori cluster 
(IHernandez et al.l 120071 ). Second, in contrast to the IC 348 population, the Cepheus region 
does not produce a rich TD population in the interval —2.6 < < —1.8. The X-ray selected 
Cepheus population has 20/139 (14%) TDs among disk-bearing stars compared to 70/163 
(43%) in IC 348. This apparent difference in TD frequency is discussed in § 16.31 



3.3. IRAC Color-Color Diagrams 

Comparison of Spitzer-IRKC source positions in the [3.6] — [4.5] vs. [5.8] — [8 . 0] dia gram 



with the result s of mo dels of stars with dusty disks and envelopes led I Allen et al.l (120041 ) and 



Megeath et al.l (20041) to propose this diagram as an excellent tool for PMS classification. 



Hartmann et al.l (l2005l ) show that similarly successful PMS classification schemes can be 
obtained at shorter wavelengths using the [3.6] — [4.5] vs. [4.5] — [5.8] and Kg — [3.6] vs. 
[3.6] — [4.5] diagrams. In Figure [Tj, we compare the SED-based classifications of Cepheus 
X-ray stars derived in §3. II with the expected loci of PMS stars in IRAC color-color diagrams. 
Here we also begin comparison of the disk characteristics in three spatial subregions: the 
unobscured Cep 0B3b cluster, the optically-bright S 155 nebula, and the Cep B molecular 
cloud core. The regions are measured radially from the location of the hot molecular core 
at a = 22^^57"^15!2, 5 = +62°37'll'.'l (J2000). 

Figure [7] shows the three IRAC color-color diagrams for the sources in the three sub- 
regions. Seventy-seven Cep 0B3b sources are missing from the plots due to inadequate 
photometry, often due to location outside of the full four-band coverage. Sources are also 
omitted from the S 155 and Cep B areas when photometry at the longer 5.8 fim and 8.0 /im 
bands are absent or inaccurate due to bright nebular emission. The positions of stars in the 
color-color diagrams generally agree with our previous classification of stars as disk-bearing 
and diskless PMS stars based on comparison with IC 348 stars ( §3.ip . shown as blue and red 
symbols in Figure [71 The F2 flag of Table [2] indicates sources with their disk classification 
from 13 . 1 1 inconsistent with the simple color criterion that diskless stars have [3.6] — [4.5] < 0.2, 
[4.5] — [5.8] < 0.2 or [5.8] — [8.0] < 0.2. Several sources with such "discrepant" locations on 
color-color diagrams may represent additional cases of transitional disks. 
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4. NON-CHANDRA INFRARED-EXCESS MEMBERS 

Due to the relatively short Chandra exposure and large distance to the Cepheus region, 
the X-ray PMS sample is only complete down to ~ 0.5 for Cep 0B3b/S 155 Class III 
stars, less complete for Cep 0B3b/S 155 Class II stars, and even less complete for Cep B 
embedded objects As the stellar Initial Mass Function peaks around ~ 0.3 Mq, hun- 

dreds of Cepheus stars should be present in the Chandra field of view with their X-ray 
luminos ities below the sensi tivity limits of the Chandra observation of logL^ ~ 29 — 29.5 



erg s~^ (IGetman et al.ll2006l ). With an age around ~ 2 — 3 Myr and younger (Appendix A), 
a large fraction of the Cepheus stars with masses below 0.5 Mq in our field will have op- 
tically thick disks detectable with IRAC. Selection by IRAC IR-excess should even detect 
many brown dwarfs in the region. Shallow IRAC exposures efficiently detected disks around 
brown dwa rfs in the Orion clou d at a distance of ~ 400 pc over several magnitudes in 



brightness (iLuhman et al.l l2008bl . and references therein). The sensitivity at the Cepheus 
distance around ~ 725 pc will be reduced by ~ 1.3 mag, sufficient to detect a considerable 
fraction of the disk-bearing brown dwarfs. Thus some fraction of non-X-ray PMS stars with 
lower- masses are expected to be detected in the IRAC survey within the Chandra field of 
view. 

Figure [Hb — c shows IRAC color-color diagrams for 774 IRAC sources in the C handra 
ACIS-I field of view that are not associated with X-ray sources. We restrict these diagrams to 
sources with photometric errors < 0.1 mag in all four IRAC bands. Using the IR color-color 
classification criteria presented in § 13. 3[ we identify 224 non- X-ray IR-excess stars. Table [3] 
gives IR properties for the 224 non-X-ray IR-excess sources. 

Figure — f compares X-ray-selected and IR-excess-selected samples in the [4.5] vs. 
[4.5] — [8.0] color-magnitude diagram. As expected, the majority of these new PMS candidates 
are fainter than the X-ray selected sample, although there is a considerable overlap in the 
distributions around 14 < [4.5] < 12 mag (~ 0.1—0.3 M©). Most of these new stars thus have 
masses in the < 0.5 M© and some in sub-stellar range. Panels show that, despite the different 
selection criteria, the color distributions of these IRAC-selected stars are similar to those of 
the X-ray selected PMS stars. However, several of the weakest and reddest ([4.5] > 14 mag 
and [4.5] — [8 . 0] > 0.5 mag) IR-selected sources may be unrelated extragalactic objects 



(IHarvey et al.l l2006l ). Four probable extragalactic contaminants from the X-ray-selected 



sample lie near this region of the diagram (grey diamonds in Figure [Hli) . 

In the unobscured Cep 0B3b subregion, the combined X-ray-selected and IR-excess- 
selected sample should give a complete census of the disk-bearing population down to 0.2 — 
0.3 M0, a nearly complete census of the disk-free population down to 0.5 M©, and a fraction 
of both the disk-bearing and disk- free populations into the brown dwarf regime (§S]). But 
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due to the limited point source sensitivity from the non-uniform IR nebular emission, the 
identified PMS populations of the S 155 and Cep B subregions are less complete. This can be 
seen from the paucity of sources fainter than [4.5] > 14 mag in S 155 and [4.5] > 13 mag in 
Cep B (Figure [S^ — /). More accurate information on mass completeness limits of Cepheus 
PMS stars is given in the next section. 



5. INITIAL MASS FUNCTION OF THE CEPHEUS CLUSTER 

In the area covered by the Chandra field and all four bands of the IRAC mosaic, we 
examine the stellar IMF separately for diskless and disk-bearing stars in each of the three 
sub-regions. For a given mass, the X-ray detection efficiency of Class III stars is somewhat 
high than that of Class II stars (Appendix [B]). To compensate for this effect, in this IMF 
and the following disk evolution analyses, we use the disk-bearing stellar sample that is the 
combination of the Chandra and non-Chandra IR-excess member young stars (§!]) while 
retaining the diskless stellar sample as purely composed of X-ray stars. In this IMF analysis 
section, we exclude X-ray stars lying outside of the 4 IRAC band mosaic, as well as the most 
massive young star in the region, 07Vn star HD 217086 {Chandra # 240) which is the only 
star with mass above 4 Mq (upper boundary of considered IMF mass range; Figure [9]). We 
also exclude 51 non-Chandra IR-excess possible members with inferred masses below 0.1, as 
some of them may still be extragalactic background objects (§!]). 

Figure [Hi shows that the mass distributions of the Cep 0B3b sub-region stars, both 
diskless (red) and disk-bearing (blue), nicely follow the shape of the Galactic field IMF, 
showing a Salpeter powerlaw slope between 0.5 — 3 M©, and a peak around 0.3 — 0.5 Mq. 
The decline seen below 0.3 — 0.5 Mq is attributable to our incompleteness limits and is 
probably not intrinsic to the cluster. The pure Chandra Class III sample is complete down 
to 0.5 Mq, while the combined Chandra and non-Chandra Class II sample is complete for 
> 0.3 Mq Cep 0B3b stars. Assuming a Galactic field IMF shape, we estimate the total 
PMS population in the Cep 0B3b sub-region to be ~ 750 stars down to 0.1 Mq. 

An importan t negative result here is that we do not confirm the tentative explanation 
by I Get man et al. I Jiooe) that the non-standard X-ray luminosity function of the Cep 0B3b 



is due to an anomalous 'bottom-heavy' IMF with an excess of 0.3 Mq stars. Although a 
slight difference between the combined Chandra and Spitzer sample and the Galactic field 
IMF is seen (Figure [Qji), it is not statistically significant. 

Figure [9^ shows that mass distributions of the S 155 sub-region stars, although with 
fewer stars than in Cep 0B3b, also follow the Galactic field shape with mass completeness 
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limits of 0.5 M© for both Class II and Class III stellar samples with an estimated total PMS 
population of ~ 200 S 155 stars down to 0.1 M©. 

The observed mass distribution of the embedded Cep B cluster may be either top-heavy 
or, more likely, simply lacking unidentified low- mass disk-free PMS stars due to Chandra 
sensitivity loss from obscuration (Figure [H:). A similar situa tion is noticed in the study of the 
embedded population of the bright-rimmed cloud IC 1396N (iGetman et al.ll2007l ). In the case 
of usual IMF with still unidentified low-mass stars, the observed Cep B stellar population 
may be complete only somewhere above 1 M©, and the total intrinsic PMS population of the 
Cep B may reach a few hundred stars. The triggered population is th us much larger than 
that identified by early radio and infrared techniques (ITesti et al.lll995l ). 



DISK EVOLUTION 



6.1. Mass Dependence 



Growing evidence has emerged fr om IRAC young; cluste r studies that more massive 
PMS stars have shorter disk lifetimes. Luhman et al. ( 2008al lbl) find that disk fractions in 
the 2 — 3 Myr old IC 348 and a Ori clusters decrease from ~ 50 — 60% in brown dwarfs to 
~ 30 - 40% in stars with M > 1 M©. In the older 5 Myr old Upper Sco and NGC 2362 
cluste rs, the disk fraction decreases from 10 — 20% in 0.1 — 1 M© stars to < 1% in higher-mass 
stars (ICarpenter et al.ll2006l : iDahm fc Hillenbrand! 120071 ) . The trend may no t be universal: 



the dis k fraction of M > 1 M© stars in the 2 Myr old Cha I cluster is > 60% (ILuhman et al. 
2008al ). much higher than that found in IC 348 and a Ori at similar ages. 



The Figure [TOl compar es the mass-dep e ndent disk fractions of Cep 0B3b and the other 
three regions examined by iLuhman et al.l (j2008bl ). Our sample here is the combined 215 
Chandra (classified as "DSK" or "NoD" in Table [2]) and 142 disk-bearing non-Chandra 
(Table [3]) Cep 0B3b sources within the overlap area of 4 IRAC band and Chandra fields, 
omitting the most massive young star in the region, 07Vn star [Chandra # 240) with 
Mj < mag, and 37 non-Chandra IR-excess possible members of Cep 0B3b with inferred 
masses below 0.1 M©. Appendix [B] shows that the relatively short Chandra observation 
is only sensitive to stars down to 0.2 M©, so we have no Cep 0B3b coverage for very 
low mass Class III stars and brown dwarfs. Cep 0B3b stars have average ages around 
2 — 3 Myr old (Appendix |A]) similar to the ages in Cha I, IC 348, and a Ori; possibly 
younger stars in the S 155 and Cep B region are omitted here. Following Luhman et al., 
the absolute J-band magnitude Mj serves as a proxy of mass. Our Mj values are obtained 
using individual Aj extinctions from Tables [2}l3], which in turn were derived through the 
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evolutionary model-dependent procedure assuming the distance of 725 pc and the age of 
3 Myr. Individual magnitude accuracies should be sufficient to reliably assign individual 
sources to the coarse, 2 mag- wide Mj bins. Erro rs on disk fraction s have been estimated using 



binomial distribution statistics as described by iBurgasser et al.l (120031 ) . For the Cep 0B3b 
stars, Mj bins have been systematically shifted by -1-0.3 mag from those of other regions in 
order to include in the second Mj bin the ~ 0.5 Mq mass stars, stars for which the combined 
Chandra and non-Chandra Cep 0B3b sample is complete (^. 

Figure fTOl shows that the disk fraction of the Cep 0B3b stars (circles) for the two highest 
mass bins is around 45% and agrees with the fractions seen in IC 348, but appears higher than 
the comparison clusters for the lowest < 0.5 mass bin. The most plausible explanation 
for the reduced disk fraction in the < 0.5 Mq bin is incompleteness of the Cep 0B3b sample 
at low masses (^. Here, our short Chandra observation is not sensitive to many low-mass 
Class III stars and the observational bias towards non-Chandra disk-bearing stars affects 
the sample selection. 

We thus find no significant dependence of disk fraction on mass in the Cep 0B3b 
Chandra plus non-Chandra selected sample above 0.5 MfJ^. This agrees with the lack of 
mass dependence seen in IC 348 and a Ori in the > 0.2 regime. The increased fractions 
reported in these clusters were found in the < 0.2 regime which is not covered by 
our Cep 0B3b selected sample. However, we do not support studies in other clusters that 
report a rapid loss of IRAC-band disk emission in intermediate-mass stars. These studies 
might consider possible sample incompleteness effects, such as an under- sampling of diskfree 
Class HI low-mass members. Our inclusion of the X-ray selected sample is specifically 
designed to be nearly free of this bias. 



6.2. Spatial Gradients 

We now consider spatial distribution of the disks around the mass-complete Chandra 
plus non-Chandra PMS samples over the area covered by the 4 IRAC mosaics and the 
Chandra field. Figure [TT] shows the spatial distribution, and Table |4] provides quantitative 
details for two mass strata, two azimuthal zones (NE and SW), three radial subregions 
(Cep 0B3b, S 155, and Cep B), and three radial layers discussed below (inner, intermediate, 
and outer). The azimuthal zones are separated by the line that roughly bisects the optically- 
bright (Figure [1]) and IR-bright nebular emission. Radial distances are measured from the 



^The results are similar with a choice of longer distance of ^ 850 pc and thus younger age of ~ 2 Myr 
for the Cep 0B3b cluster (Appendix IXj) . 
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hot core of the molecular cloud (a = 22^'57™15!2, 5 = +62°3ril'.'l (J2000)). Below we 
consider the Cepheus sample of 220 Chandra plus non-Chandra PMS stars with stellar 
masses above 0.5 Mq, the mass completeness limit for Cep 0B3b/S 155 stars as well 
as its mass-stratified sub-samples. The observed stellar population of the embedded Cep B 
cluster may be complete somewhere above 1 M0 and thus the analysis of the > 1 Mq 
stellar sub-sample may provide more realistic disk fraction estimates for this sub-region. 

Examining first azimuthal dependence, we find that the NE zone of the region has 
~ 1.5 — 1.7 times higher stellar surface density than the SW zone. There is also an indication 
for possible slight increase in disk fraction from NE to SW zones of the Cep 0B3b/S 155 
region, although the effect is not statistically significant (Table H]). 

Examining next the radial dependence, a dramatic trend is seen: the disk fraction 
increases from 40 — 50% in the Cep 0B3b subregion to 50 — 60% in the S 155 subregion 
to 70 — 80% in Cep B (Table H]). This trend is elucidated in more detail in Figure [12] as a 
running average disk fraction. Here the running disk fraction for the > 0.5 Mq mass stratum 
(red) is evaluated within a 2' wide sliding window at 1' intervals, while the 0.5 < M < 1 Mq 
(black) and M > 1 Mq (blue) mass strata are evaluated within a 4' wide sliding window at 
2' intervals. Note these values are not independent because a given star is included in two 
adjacent values. The effect could either be viewed as a smooth trend in disk fractions, or 
a step-function at 30% (outer layer), 60% (intermediate layer) and 70 — 80% (inner layer) 
disk fractions in 0.6 — 0.8 pc radial bins (grey bars). It is important to note that the outer 
edge of the second bin (intermediate layer) is coincident with the projected position of the 
ionizing source of the region, 07V star HD 217086. 



If one adopt the simple relation between stellar age and disk fraction derived by lHernandez et al 



(120071 . see their Figure 14), the Cepheus disk fraction gradient corresponds to an age gradient 
from 3 — 5 Myr in the outer layer (part of Cep 0B3b subregion), 2 — 3 Myr in the inter- 
mediate layer (part of Cep 0B3b/S 155 region), and ~ 1 Myr in the inner layer (embedded 
Cep B subregion). This inferred age gradient agrees with that estimated from stellar {vs. 
disk) properties in Appendix Rl 



6.3. Transition Disks 



As noted earlier in § 13.21 the Cepheus 0B3b population appears deficient in transition 
disks compared to the 2—3 Myr old cluster IC 348. We in vestigate this further here, including 
comparison to the low-r aass 1 Myr old Coronet cl uster (jSicilia-Aguilar et al.ll2008l ) and the 
4 Myr old Tr 37 cluster (jSicilia-Aguilar et al.ll2006l ). Tr 37, like our Cep 0B3b field, harbors 
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an 07-type star. 



It is important to recall that different studies of different star forming regions use dif- 
ferent criteria for identifying TDs. In the Tr 37 and the Coronet cluster studies, TDs are 
identified with a photometric excess longward of ~ 6/im. In the stud y of IC 348, small 
excesses around 2 /xm define "anemic disks" systems (ILada et al.l 120061 ) . Our classification 
of Cepheus stars here relies on an intermediate criterion of excess longward of ~ 3. 6/im. 
Thus, if the fraction of transition disks were intrinsically the same in all clusters, we expect 
the reported fractions to be largest in IC 348, intermediate in our sample, smallest in Tr 37 
and the Coronet. Furthermore, inter-cluster comparison of TD fractions should also con- 
sider similar mass strata, as lower-mass PMS sta rs may have flatter disks with more frequent 



inner holes than more massive PMS star s (e.g. iHartmann et al.l l2006l : ISicilia-Aguilar et al. 



20081 ) . From Table 2 of iLada et al.l (120061 ) . a TD fraction among disk-bearing stars in IC 348 
decreases from 0.45 ± 0.04 (64/142) for M-type stars to 0.3llo:o9 (4/13) for G-K-type stars, 
although this decrease is not statistically significant. The Tr 37 study combined stars from 
late G to M2-type (~ 0.4 - 2 M©) and reported a 0.10l[5:[J^ (11-14 out of 140 - 150) TD 
fraction among disk-bearing stars. The Coronet cluster study included only M0-M8 objects 
and reported a 0.50+°;^^ (7/14) TD fraction. 

Our classification of X-ray selected PMS stars in the Cep B/Cep 0B3b obtained 20 
TDs out of 139 disk-bearing stars giving a TD fraction of 0.14]'^q'q2- Most of these (13 out of 
20) are in the unobscured Cep 0B3b subregion and none are present close to the Cep B hot 
core. To perform a mass-stratified analysis, we consider here only the most rich star sample, 
Cep 0B3b. From Table [21 we find the Cep 0B3b TD fraction may be mass-dependent with 
0.29lo:o7 (8 out of 28) among M < 0.5 Mq (M-type) disk-bearing stars and 0.1llo;o3 (5 out 
of 47) around more massive stars. 

These results are roughly consistent with the previous studies, once the TD definitions 
and mass ranges are taken into account. Our overall ~ 14% TD fraction is intermediate 
between the ~ 40% reported for IC 348 and ~ 10% reported for Tr 37, as expected for an 
intrinsically constant TD fraction. Our finding that M stars have a higher TD fraction than 
more massive stars agrees with the indication for trend seen in IC 348 and the high fraction 
found in Coronet M stars. We thus conclude that the fraction of Cep 0B3b TDs agrees with 
that of other clusters despite the apparent deficit noted earlier in Figure [61 
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7. IMPLICATIONS FOR TRIGGERED STAR FORMATION 

7.1. Analogy with Bright Rimmed Clouds 

The observed picture of the Cep B region is in many respects reminiscent, on a larger 
scale of several parsecs, of those of smaller bright-rimm ed clouds (BRCs) an d cometary 
globules (CGs) found on the edges of giant H 11 regions (ISugitani et al.lll99ll ). BRCs are 
isolated clouds surrounded by ionized rims facing the exciting star(s) with their dense cores 
close to the rims. BRCs are modeled as externally illuminated, photoevapor ated and ablate d 
into elongated head-tail morphologies by ultraviolet radiation of OB stars (lReipurthlll983l ). 
It is likely that pressure from the ionization shock front at the surface propagates through a 
globule and overcomes the magnetic, turbulent and thermal pressure that supports it against 
collapse, thereby triggering localized star formation. This astrophysics of r adiation-driven 



implosion (RDI) in molecular globules h as been extensively studied (e.g., iBertoldil Il989 
Lefioch &: LazarefBll994j : iMiao et al.ll2009l ). RDI models typically consider the triggering of 
a single star formation episode in a small (< 1 pc) cloud with tens of Mq of molecular gas. 
The characteristic timescale for producing cometary morphologies and inducing gravitational 
collapse varies with initial conditions from 0.1 to ~ 1 Myr. 

In mo st cases, molecular, IR, B.a surveys of BRCs trace only the most recently formed 



stars (e.g., ISugitani et al.lll995l : lOgura et al.ll2002l : [Thompson et al.ll2004bl : lUrquhart et al. 



20061). In a few cases, Chandra observations have added the diskfree PMS popu lations 
(ICetman etaPboO?! : Isanchawala etaPboO?! : Icetman et al.ll2008al : IWang et al.ll2009 h. Trig- 
gered BRCs often show a few embedded mid-IR sources denoting protostars, while Ha, JHK 
and X-ray surveys reveal small clusters of disk-bearing PMS stars within and in front of the 
bright rim. In a few cases, spatial-age gradients in the stellar population are seen where 
the youngest stars are embedded and older stars are aligned tow ard the ionizing sources 
( iMatsuyanagi et al.l l2006l : lOgura et al.l 120071 : ICetman et al.l 120071 ). This directly supports 
the RDI mechanism and implies that the existing clouds have been actively forming stars 
for several million years. 

The Cepheus region possesses all of the observational features of an RDI triggered star 
formation region: 



1. The presence of exciting star(s) and a mole cular cloud surrounded by an ionized rim 
facing the exciting stars shown in Figures[Tll2]( Minchin et al1ll992 : Beuther et al. 2000 ): 



2. The p resence of a dense molecular core close to the rim (lYu et al.lll996l : iBeuther et al, 
2000h : 
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3. The spatio-temporal gradient of young stars oriented towards the exciting star(s) shown 
in Figures 

The major difference of the Cepheus region from other BRCs is its larger scale: its extent is 
~ 4 pc linear size, its molecular mass is ~ 200 Mq, its observed PMS X-ray stellar population 
is > 60 embedded and up to ~ 300 older unobscured members, as well as estimated intrinsic 
population of ~ 1000 unobscured and possibly a few hundred embedded stars (^. The large- 
scale morphology is also different from most other BRCs. Rather than lying on the edge of a 
circular HII region surrounding a concentrated OB association, it is an extension of a giant 
molecular cloud protruding into a large evacuated region produced by a partially dispersed 
OB association. Another unusual case of a very large BRC with st ar formation triggere d 



over several million years is the high-latitute cometary globule CG 12 (IGetman et al.ll2008al ). 



7.2. Implications for the Radiative Driven Implosion Mechanism 

Based on the RDl triggered star formation concept, we can propose some new in- 
sights into the star formation process of the Cepheus region and make some simple model- 
independent estimates for a shock propagation velocity and a star formation efficiency (SFE) 
of the RDI process. In §6.21 and Appendix A, we establish the spatio-temporal gradient of 
Cepheus PMS stars from the center of the cloud towards the HD 217086 07Vn star. Assum- 
ing a distance around 725 pc, the innermost 0.6 — 0.8 parsec adjoining the hot core of the 
cloud has the youngest (< 1 — 2 Myr) stars with the highest seen disk fraction (70 — 80%). 
The intermediate 0.6 — 1 parsec layer around S 155 has older (2 — 3 Myr) stars with an 
intermediate disk fraction (~ 60%) (Figure [12]). The outer layer of mostly diskless stars in 
the unobscured Cep 0B3b cluster has the lowest disk fraction (~ 30%)|2| 

We consider two scenarios for RDl-induced star formation in this region. First, a rel- 
atively slow shock with speed < 0.5 km s~^ may have passed through the molecular cloud, 
triggering star formation in the outer layers around 2 — 3 Myr ago and continuing to trigger 
star formation close to the Cep B hot core today. This slow shock speed is consistent with the 
0.6 km s~^ propagation rate inferred from the age gradient of the PMS stars in the globule 



IC 1396N, located ~ 11 pc projected distance of its ionizing 06e type star (IGetman et al. 



20071 ). In this scenario, stellar kinematic drift from their birthplaces plays an important role 



in relocating PMS stars. The 07Vn star HD 217086 star is probably the principal ionization 



■^If the distance to the cloud is ~ 850 pc, the ages of the Cep B and S 155/Cep 0B3b subregions are older 
by 0.5 — 1 Myr but the 1 Myr difference between them remains (Appendix A). 
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source but other OB stars (Figured]), including the BlVn star HD 217061 currently located 
within the SW zone of the S 155 subregion, probably also play a role. If the original cloud 
were only slightly larger than the size observed today, the RDI shock would have slowly 
propagated through ~ 1 pc of the Cep B molecular cloud over ~ 2 — 3 Myr years, implying 
a shock velocity around < 0.5 km s~^. The stars formed early in this process would have 
drifted 2 — 10' (0.4 — 2 pc) in all directions, populating much of the Chandra field of view 
with triggered stars. Such star drifting would require a stellar velocity dispersion around 
1 km s ~^. This is expected from th e turbulent velocities within a cloud several parsecs in 
extent (lEfremov &: Elmegreenlll998l ). whic h can be supplemented by dyn amical interactions 
between protostars during star formation (jBatdl2008l : iFiiresz et al.ll2008l ). 

A second scenario involves the passage of a faster shock propagating at ~ 1 km s~^ 
through a cloud that was originally much larger than the Cep B cloud we see today. In 
this model, most of the original cloud material has ablated and most of the stars in the 
Chandra field formed in the cloud during the past ~ 2 — 3 Myr. Present locations of the 
intermediate layer around S 155 and the inner layer around the molecular core (separated by 
0.8 — 1 pc. Figure fT2l) may be directly associated with the passage of the RDI shock ~ 2 — 3 
and < 1 — 2 Myr years ago, respectively. Star drift plays a less important role here. A similar 
fast shock speed is observed in the triggered star forn aation in the molecular pillars of the 



Eagle Nebula (M 16) and other bright rimmed clouds (IFukuda et al.l 120021: [Thompson et al. 



2004bl ). and is consistent with the theoretical models of iMotoyama et al.l ( 120071 ). 



A critical discriminant between these scenarios is the SFE of the cloud. Integrating the 
stellar masses in Figure [H the total stellar mass of the Cep 0B3b disk-bearing population 
assuming a standard Galactic field IMF is ~ 130 M© in our observed field. With the average 
disk fraction of 45% ( §6.H Table Hj), we add ~ 160 M© of diskfree stars and ~ 20 Mq from 
HD 217086 to give a total mass of the Cep 0B3b population around 310 Mq. A similar 
estimate for the S 155 population from information in Figure [HI and Table H] gives a total 
stellar mass around 80 Mq. Assuming that the observed Cep B embedded population is 
complete only somewhere above 1 M© (^ and follows a standard Galactic field IMF as 
well, we use star count and disk fraction information for the highest mass bin from Figure [9] 
and Table m to obtain an estimated total mass of the Cep B population of 120 — 250 Mq. 
The mass of the observed molecul ar gas in the Cep B molecular por tion of the Cepheus 
molecular cloud is only ~ 200 Mq JYu et al.l 1 19961 : iBeuther et aPboOoh . 

The inferred SFE for the embedded Cep B population alone is thus around 35% — 55%, 
at the top of typical SE E range of 10 — 40% measured in other active star forming regions 
(lElmegreen et al.l |2000| . and references therein) . This is far above th e SFE found in smaller 



molecular globules triggered by UV shocks (e.g. iGetman et al.ll2007l ). For the total Cepheus 
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stellar population within the Chandra field, the apparent SFE is extremely high around 
> 70%. This is unrealistically high and indicates that most of the Cep 0B3b stars were not 
formed from the presently seen cloud material but arose from an earlier generation of star 
formation from gas that is no longer present. This supports the second scenario involving a 
fast shock passing through a much larger and more massive original cloud. 



Two other aspects of RDI triggered star formation can be discussed. First. ISugitani et al. 



(119911 ) provided observational evidence for a non-standard IMF where intermediate-mass 
stars are preferentially formed over lower mass stars in BRCs. In the case of the embedded 
triggered PMS population of IC 1396N, iGetman et al.l (120071 ) also found indications for a 
nonstandard IMF biased towards higher mass stars, but cautioned about possible observa- 
tion selection effects. A similar situation arises here in the case of the embedded Cep B 
population While the shape of the Cep B stellar mass distribution is not definite yet, 
the stellar mass distribution of likely RDI triggered stellar population of the Cep 0B3b/S 155 
region agrees well with a standard Galactic IMF, at least down to ~ 0.2 — 0.5 M© (^. 

Second, most theoretical calculations of RDI triggering involve small globules and a 
single episode of triggering. The 2—3 Myr range of st ellar ages found in the Cep B/Cep 0B3b 
region, and an even wider age spread found in CG 12 (IGetman et al.ll2008al ). implies repeated 
or continuing star formation over millions of years when the RDI mechanism occurs in a larger 
molecular cloud. 



8. CONCLUSIONS 

We present a Spitzer IRAC observation of the Cepheus B molecular cloud, the S 155 HII 
region on its periphery, and a portion of the nearby Cep 0B3b OB association. The goals of 
this work are to study disk evolution of the nearly disk-unbiased combined Chandr a- ACIS 
and Spitzer-IKAC selected samples of PMS stars and to provide new clues of the sequential 
triggered star formation process in the region. 

Out of ~ 400 X-ray emitting PMS stars in the region, 354 are classified as diskless or 
disk-bearing stars based on IR photometry: 161, 36, and 18 (75, 26, and 38) diskless (disk- 
bearing) stars in the Cep 0B3b, S 155, and Cep B subregions, respectively. We immediately 
see that samples selecting only IR-excess stars miss the majority of stars outside the Cep B 
molecular core. Among all X-ray emitting disk-bearing systems, only 4 are Class I protostars; 
they lie in the younger embedded Cep B cluster ([|3]). In addition to the ~ 400 X-ray emitting 
PMS stars in the region, we identify > 200 non-Chandra IR-excess low-mass members of 
the region (§!]). For the lightly-obscured Cep 0B3b and S 155 clusters around the cloud, the 
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combined Chandra and non-Chandra PMS samples are complete down to 0.5 M©, while for 
the embedded Cep B cluster the observed PMS sample may be complete only above 1 Mq 

The spatial distribution of the disks around stars of the combined Chandra — Spitzer 
PMS samples reveals a dramatic spatio-temporal gradient: younger stars are clustered 
around the hot core of the molecular cloud, while older stars are dispersed in the direc- 
tion of the primary ionizing 07-type star HD 217086 ( §6.21) . Disk fractions fall from 80% 
to 30% across the S 155 HII region. This finding constitutes strong endorsement for the 
basic model of triggered star formation in the region by passage of shocks driven by the 
OB stellar ionization and/or winds into the Cep B molecular cloud. Our estimate of the 
shock velocity of ~ 1 km s~^ agrees with empirical estimates obtained for a number of other 
shocked bright-rimmed clouds and with theoretical models. The IMF of the Cep 0B3b/S 155 
triggered stellar population exhibits a standard Galactic field shape; we do not confirm our 
earlier suggestion of a bottom-heavy IMF. The unusually large apparent value of star forma- 
tion efficiency of > 70% based on the present day molecular cloud indicates that most of the 
Cep B molecular material has likely already been evaporated or transformed to stars. Con- 
trary to the current RDI simulations of typical small globules which model a single triggering 
event over < 1 Myr time period, our observational findings support a picture of multiple 
episodes of triggering over millions of years (©. 

Using the nearly disk-unbiased combined Chandr a- ACIS and Spitzer-IRAC selected 
stellar sample, we also obtain several other results. First, the disk fraction as a function of 
stellar mass in the Cep 0B3b cluster is similar to that of the IC 348 PMS stars with ages 
similar to Cep 0B3b ( §6.ip . Second, the fraction of Cep 0B3b transition disks agrees with 
that of other clusters, such as Tr 37, Coronet, and IC 348 ( §6.3p . Third, in agreement with 
the disk fraction analysis of §6.21 we obtain photometric ages of < 1 — 2 and 2 — 3 Myr for 
the embedded Cep B and lightly-obscured Cep 0B3b clusters, respectively using the optical 
color-magnitude diagram (Appendix and assuming a distance of 725 pc. Fourth, the 
effect of reduced X-ray emission in disk-bearing versus diskless stars which was originally 
established in studies of ONC and Taurus regions is confirmed in the Cepheus systems 
(Appendix [B]) . 



A. APPENDIX A: DISTANCE AND AGE OF THE CEP B/CEP OB3 

COMPLEX 



The distance to the Cep B/Cep 0B3 complex is somewhat un certain. From the op tical 
photometry of the early-type members of the Cep 0B3 association, iBlaauw et al.l (Il959l ) and 
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Craw: 



ord &: Barnea (Il970l ) obtained a distance estimate of 725 pc. Later lMoreno-Corral et al 



( 119931 ) combined near-IR photometry with the earher optical photometry and spectroscopy 
data of bright members to derived a distance estimate of 850 pc. A recent VLBI parallax 
measureme nt using the methano l maser in the Cep A molecular core gives a distance of 
~ 700 pc ( Moscadelli et all 2009). A distance of 870 pc has been infe rred for the nearby 
Cep 0B2 association (IContreras et al.l l2002l : ISicilia-Aguilar et al.ll2005l ). Our X-ray study 



(I Get man et al 



by 



Pozzo et al. 



20061) adop ted the 725 pc estimate while the optical studies of the Cep 0B3b 



(120031 ) and lMayne et al.l (120071 ) adopted a distance of 850 pc. 



The chosen distance influences age estimates of the lower-mass stars as PMS evolution- 
ary isochrones on the Hertzsprung-Russell diagram (HRD) depend on absolute luminosities. 
Early photometric age estimates of the OB population indicated a relatively old cluster 

For the X-r ay PMS stars in the 



199 



around ~ 4 — 5.5 Myr (iBlaauwl 119641 : iJordi et al 
Chandra field, we adopted an age of ~ 1 M yr (iGetman et al 



2006) based on the HRD 



of five PMS X-ray sources (jPozzo et al.ll2003l ). Recently, iMayne et al.l (120071 ) obtained V- 
and Jc-band photometry data for an area including most of the Chandra ACIS-I field and 
measured a ~ 3 Myr photometric age for the low-mass Cep 0B3b population. It is unclear 
from any of these earlier datasets whether the age is uniform across the region. 

We seek here to estimate stellar ages from the V- and /r-band colo r -mag nitude diagram 
for the Chandra X-ray PMS stars using the photometry of lMayne et al.l (120071 ) listed in Table 
H The 2MASS n ear-IR photometry was first used with the comb ination of the synthesized 



Siess et al.l J2000h (for 1.4 ^ M ^ 7.0 Mq) and lBaraffe et all Jl998h (for 0.02 ^ M ^ 1.4 Mq) 



PMS evolutionary isochrones to obtain individual extinction estimates for Chandra sources. 
We obtained different sets of source extinctions {Av^nir) us ing the J versus J — H color- 
magnitude diagram (similar to the one shown in Figure 5 of iGetman et al.ll2006l ) with two 
trial assumed distances (725 and 850 pc) and three trial assumed ages (1, 2, and 3 Myr). 
The six sets of source extinctions were applied to correct the observed optical magnitudes 
of Chandra sources and to place them on the V versus V — Ic diagram. Based on the 
val ues of E{B — V) = 0.91 mag and E{V — I) = 1.18 mag at Ay = 2.8 mag adopted 
by iPozzo et al. J2003h and the reddening laws of IWinklerl (Il997l ). we used the reddening 
relationships Aj^/Ay = 0.58, Aj/Ay = 0.27, and Aff/Ay = 0.16. The Cep B optical stars 
have an average extinction < Ay^^jji >~ 7 mag while S 155 stars have < Ay^jyui >~ 3 mag 
and Cep 0B3b stars have < Ay^^iR >~ 2 mag. Within each of the areas, the extinctions of 
disk-bearing stars are systematically higher than those of diskless stars. 

As the inferred visual source extinctions differ by Ay < 0.5 — 1 mag for the six trial 
combinations of distance and age, and the reddening vector on the optical color-magnitude 
diagram is almost parallel to the PMS evolutionary tracks, the X-ray sources occupy very 
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similar loci in the V — Ic color-magnitude diagram for all six cases. The X-ray sources 
are centered mostly around the ~ 2 — 3 Myr isochrons assuming the distance of 725 pc, 
or the ~ 2 Myr isochrone assuming the distance of 850 pc. Figure [12] shows the color- 
magnitude diagram for 725 pc. The plotted symbols stratify the 286 X-ray sources with 
optical photometry according to location in the star forming region and the presence of 
disks: 62 disk-bearing (143 diskless), 22 disk-bearing (31 diskless), and 17 disk-bearing (11 
diskless) for the Cep 0B3b, S 155, and Cep B subregions, respectively. 

In Figure[13l about thirty X-ray sources have V — {V ~Ic) locations inconsistent with any 
PMS model tracks. These outliers are not foreground or background stellar contaminants, 
as more than half of them have optically thick disks and very few field stars are predicted 
to be captured in the Chandra image. Most of these outliers have unusually high near- 
IR extinctions Ay^Mm among the X-ray young stars seen in the optical bands. For the 
majority of the outliers, the Av^nir values are comparable (within 30%) to the absorption 
Av,xray inferred from the X-ray median energy, MedE. MedE is first converted to equivalent 
hydrogen column densit y assuming solar abu ndances using the MedE — Nh relation found 



for the COUP sample (iFeigelson et al.l 120051). and N m is converted to Ay^xray assuming a 
gas-to-dust ratio of ~ 2 x 10^^ cm~^/mag (jRyter 1996 ). This may point to discrepant optical 



magnitudes and/or unusual conditions in these PMS systems; for example, some may posses 
disk at high inclination where the optical light is enhanced by scattering above the disk. 
Alternatively, these outliers m ay have ages ove restimated by incorrect consideration of birth 



line and accretion effects (e.g. iHartmamj |2003| ) 



From a compari son of soft X-ray absorptions and fractions of i^^s-band excess sources. 



Getman et al.l (120061 ) presented evidence that the embedded Cep B cluster is younger than 
the Cep 0B3b population within the Chandra ACIS-I field. Figure [13] supports this age 
difference. Ignoring outliers mentioned above, about half of the Cep B sources (blue symbols) 
are located closer to the 1 Myr track compared to only about one-sixth of the Cep 0B3b (red 
symbols). However, the Cep B sample is small because most of the X-ray selected members 
are embedded and do not have optical counterparts. The result supports the youthfulness 
of the Cep B population seen in the disk fraction discussed in § 16.21 

We thus cannot arrive at a clear conclusion on the distance to the star forming complex. 
However we do establish a dependence of the inferred average age on the assumed distance. 
In this paper, we adopt (not derive) a distance to the Cep B/Cep 0B3b complex of 725 pc, 
supported by the recent VLBI parallax measurement of Cep A. We therefore adopt an age 
of 2 — 3 Myr for the Cep 0B3b stars and < 1 — 2 Myr for the embedded Cep B stars. If the 
distance is closer to 850 pc, the inferred ages would be decreased by ~ 0.5 — 1 Myr. 
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B. APPENDIX B: X-RAY ACTIVITY, STELLAR MASS AND DISKS 



Originally found by lFeigelson et al.l (119931 ) from ROSAT data, it is now empirically well- 
established that PMS X-ray luminosities are strongly correlated with stellar mass, volume 
and surface area. The oc M^'^ relationship extends over 3 orders of magnitude range in 
X-ray luminosity. The clearest relationships are se en in the Chandra Orion Ultradeep Projec t 



(COUP) observation of the Orion Nebula Cluster (jCetman et al.ll2005l: iPreibisch et al.ll2005r) 



and i n the XMM-Newton Extended Survey of Taurus fXEST. lGiidel et al.lboOTi : 



Telleschi et al, 



20071 ). The astrophysical cause of this relationship is poorly understood, but probably is due 
to saturation of the magnetic dynamo in the fully convective stellar interior or on the surface 
of PMS stars. 

Figure [H] shows the — M relationship for 222 unobscured Cep 0B 3b PMS stars 



compa red to 457 lightly-absorbed [Ay < 5 mag) Orion COUP stars from iGetman et al. 



( 120051 ). Here, represents the quantity L^ c, the X-ray luminosity computed in the 0.5 — 
8 keV band and corrected for obscuration. A similar, roughly linear, relationship is seen 
in the Cep 0B3b sample, but somewhat flatter than the relationship present in the COUP 
sample. This can be attri buted to truncation effects: the Cep 0B3b observation is complete 
to logL^ ~ 29.5 erg s~^ (IGetman et al.ll2006l ). while the COUP observation is complete to 
log ~ 28.0 erg s~^. Examination of Figure shows that the Cep 0B3b X-ray sample 
alone is nearly complete down to 0.5 Mq and ~ 1 Mq for the Class III and Class II stars, 
respectively (see Figure [9]). Much of the scatter in both relations can be attributed to errors 
in COUP and Cep 0B3 stellar masses; less scatter is seen in the XEST — M plot where 
masses are carefully evaluated from accurate spectroscopy. X-ray members of the embedded 
Cep B cluster show a similar — M relationship with the completeness limit above 1 Mq. 

A second relationship between X-ray emissio n and accretion was o riginally found from 
Chandra HRC (High Resolution Camera) data (IFlaccomio et al.l l2003l ) and has been con- 
firmed in the COUP and XEST samples. This is thought to arise from a suppression of 
time-integrated X-ray emission in accreting versus non-accreting PMS systems. This is a 
more subtle effect than the X-ray dependence on mass with an amplitude of only a factor of 



2 within the 10^ range in PMS X-ray luminosities (IPreibisch et al.ll2005l : iTelleschi et al. 



20071 ) . The astrophysical cause of the mild suppression of X-ray emission in accreting PMS 
systems is also uncertain. Preibisch et al. suggest that X-r ay emission cannot a rise in mag- 
netic field lines that are mass loaded with disk material. Jardine et al. J200d)_ argue that 



the ou ter magnetosphere of accretors is stripped by interaction with the disk. iGregory et al. 



(I2007D propose that so ft X-ray emission is attenuated by dense material in accretion columns. 
Getman et al.l (l2008d ) find indications of shorter X-ray flare durations in Class II systems 
that may be due to distortion and destabilization of magnetic loop structures in accreting 
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systems. 

This effect is shown in Figure [T3] with running medians of accreting and non-accreting 
stars for the lightly-obscured Cep 0B3b and COUP samples. The running medians are cal- 
culated within a sliding window with width 0.3 de x in log(M/M(T)). Th e errors on medians 



are median absolute deviations (see Appendix B of lGetman et al.ll2008bl ). The sample selec- 
tions are not identical due to different X-ray sensitivity limits and due to different limited 
spectroscopic/photometric measurements in Orion and Cep 0B3b regions. We show with 
the • symbol 62 Cep 0B3b stars classified as disk-bearing 'DSK' systems in § [31 and 142 
COUP stars classified as active accretors having Ca II 8542 A line in emission with equiv- 
alent width of EW{Cei II) < —lA. The x symbols indicate 160 Cep 0B3b stars classified 
as diskless 'NoD' and 315 weakly accreting or non-accreting COUP stars with absorption 
equivalent width of EW{Ca II) > lA. Thirteen Cep 0B3b transition disk stars are marked 
by squares. 

The factor of ~ 2 offset between the accreting and non-accreting stars in the COUP 
samples is clearly seen to be superposed on the dominant L^ — M relationship, but the effect 
is only marginally present in the Cep 0B3b sample. The weakness of the offset here can be 
attributed to two effects. First, the incompleteness at lower-L^; values discussed above has 
a greater impact on the weaker-L^^ accreting stars than the stronger-L^; non-accreting stars. 
Second, due to the absence of a spectroscopic survey of Cep 0B3b stars, we classify stars 
here using an infrared photometric disk indicator rather than a direct accretion indicator. 
It has been established within the COUP sample that the pre sence of a disk itself d oes not 



suppress X-ray emission, but active accretion must be present (jPreibisch et al.ll2005l ). These 
two effects plausibly explain the less pronounced observed effect of the Cep 0B3b Class II 
X-ray emission suppression compared to that of the COUP stars. 
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0. 


,02 


AAAOOO 


0.2 


12 


.15 ± 0. 


,02 


11, 


,92 


±0, 


.02 


11 


.71 ± 0, 


,03 


10.99 ± 0.03 


4444 




287 


0.4 


17, 


76 


± 


0. 


,01 


2 


.40 ± 0. 


,01 


0.5 


13.38± 






12.54± ■ 




12, 


,26 


± 


0. 


,03 


UUAOOc 


0.6 


12 


.02 ± 0, 


.02 


11 


.96 


± 0, 


.03 


12 


.18 ± 0, 


18 




3349 




288 


0.0 


19, 


,09 


± 


0. 


,01 


3, 


.05 ± 0. 


,01 


0.1 


13.69 ± 


0, 


,02 


12.50 ± 0. 


,02 


12, 


,06 


± 


0. 


,02 


AAAOOO 


0.2 


11 


.68 ± 0, 


,02 


11, 


,62 


± 0, 


.02 


11 


.78 ± 0, 


,09 




4449 




289 


1.3 


12, 


55 


± 


0. 


,01 


1 


.05 ± 0. 


02 


1.6 


10.79 ± 


0, 


,02 


10.49 ± 0. 


,02 


10, 


,40 


± 


0. 


,02 


AAAOOO 


1.6 


10, 


.33 ± 0, 


,02 


10, 


,32 


± 0, 


.02 


10, 


.34 ± 0, 


,03 


10.37 ± 0.03 


4444 




290 


0.2 


19, 


,33 


± 


0. 


,01 


3, 


.03 ± 0. 


,01 


0.3 


14.63 ± 


0, 


,06 


14.03 ± 0. 


,05 


13, 


,80 


± 


0. 


,06 


AAAOOO 


0.3 


13, 


.46 ± 0, 


,04 


13, 


,30 


± 0, 


.05 










2299 


00 
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291 


0.5 


21, 


,32 


± 


0. 


,01 


3, 


.72 ± 0. 


,02 


0.5 


15.39 ± 


0, 


,07 


14.49 ± 0. 


,06 


14, 


13 


± 


0. 


,07 


AAAOOO 


0.3 


13, 


.81 ± 0, 


,02 


13. 


,69 


±0, 


.02 


13, 


.73 ± 0, 


,06 


13.17 ± 0.09 


4444 


292 


0.4 


19, 


26 


± 


0. 


,01 


2 


.99 ± 0. 


,01 


0.4 


14.06 ± 


0, 


,04 


12.93 ± 0. 


,04 


12, 


53 


± 


0. 


,04 


AAAOOO 


0.3 


12 


.16 ± 0, 


,02 


12. 


,13 


±0, 


.02 


11 


.62 ± 0, 


,05 




4449 




293 


0.4 


18, 


,02 


± 


0. 


,01 


2 


.66 ± 0. 


,01 


0.4 


13.48 ± 


0, 


,04 


12.51 ± 0. 


,03 


12, 


24 


± 


0. 


,03 


AAAOsO 


0.3 


11 


.98 ± 0, 


,02 


11. 


,94 


±0, 


.02 


11 


.88 ± 0, 


,03 


11.98 ± 0.03 


4444 




294 


0.3 


15, 


,88 


± 


0. 


,00 


2 


.01 ± 0. 


,01 


0.5 


11.71 ± 


0, 


,02 


10.93 ± 0. 


,02 


10, 


,62 


± 


0. 


,02 


AAAOOO 


0.3 


10, 


.35 ± 0, 


,02 


10, 


,34 


±0, 


.02 


10, 


.34 ± 0, 


,03 


10.27 ± 0.03 


4444 




295 
296 


0.1 


20 


11 


± 


0. 


01 


3, 


.49 ± 0. 


01 


0.2 


14.21 ± 





,04 


13.10 ± 0. 


,03 


12 


,62 


± 


0. 


,03 


AAAOOO 


0.5 


12 


.29 ± 0, 


,07 


12 


,07 


± 0, 


.07 


11 


.86 ± 0, 


,07 




2229 
9999 





Note. — This tabic is available in its entirety in the electronic edition of the journal. A portion is shown here for guidance regarding its form and content. Column 1: X-ray source number. For 
X-ray source positions see Table 1 of lGetman et al.| j2006| ) . Columns 2-4: Optical-X-ray positional offset, and optical V, — Ic magnitudes. Optical data are from fMayne et al. ( 2(]0 7j) . C olumns 
5-9: 2MASS-X-ray positional offset, 2MASS JHKg magnitudes, and 2MASS photometry quality and confusion-contamination flags. For 2MASS source names, see Table 2 of .Getman ct al.| J200§) . 
Column 10: IRAC-X-ray positional offset. For the most of the sources the reported offset is for an IRAC source from the 3.6 /zm-band image; when not available, the 4.5 /xm-band position is used. 
Columns 11-14: IRAC magnitudes derived in this work. Column 15: Four digit flag (one for each IRAC band) giving photometric apertures and level of source contamination from nearby sources 
and nebular IR emission: 2 pixel aperture with likely high level of contamination, 3 pixel aperture with likely moderate contamination, 4 pixel aperture with likely low level of source contamination. 
A '9' indicates the inability to derive photometry for one of several reasons: out of IRAC channel field of view (this is further clarified by the flag F3 from Table 2), weak or absent source below 
the detection threshold, strong contamination or confusion due to nearby source(s) or bright nebular emission. 



Table 2. Membership and Classification of X-ray Sources 



No 




cto 


Nb 


Aj 
(mag) 


M 
(Mq) 


Class 


Fi Fa 


Subregion 


F3 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 


(10) 


280 


-1 


49 ± 0.03 


4 


0.6 


0.2 


DSK 


000 


CepOB3b 


0000 


281 


-2 


60 ± 0.02 


4 


0.5 


0.5 


NoD 


000 


CepOB3b 


0000 


282 


-1 


53 ± 0.02 


4 


0.4 


0.8 


DSK 


000 


CopB 


0000 


283 


-2 


67 ± 0.03 


4 


0.5 


0.7 


NoD 


000 


CepOB3b 


0000 


284 


-2 


65 ± 0.08 


4 


0.5 


0.4 


NoD 


000 


CepOB3b 


0000 


285 


-1 


29 ± 0.04 


4 


0.6 


0.3 


DSK 


000 


CepB 


0000 


286 


-1 


73 ± 0.01 


4 


0.5 


0.5 


DSK 


000 


S155 


0000 


287 


-2 


72 ± 0.08 


3 


0.1 


0.5 


NoD 


300 


CepB 


0000 


288 


-2 


71 ± 0.06 


3 


1.0 


1.0 


NoD 


300 


S155 


0000 


289 


-2 


88 ± 0.01 


4 






FRG 


333 


CepOB3b 


0000 


290 


-2 


27 ± 0.23 


2 






FRG 


333 


CopB 


0000 


291 


-2 


41 ± 0.05 


4 


0.5 


0.2 


DSK 


TD Oil 


CepOB3b 


0000 


292 


-2 


46 ± 0.06 


3 


0.8 


0.7 


DSK 


TD 301 


S155 


0000 


293 


-2 


78 ± 0.01 


4 


0.5 


0.8 


NoD 


000 


CepOB3b 


0000 


294 


-2 


77 ± 0.01 


4 


0.7 


2.1 


NoD 


000 


CepOB3b 


0000 


295 


-2 


03 ± 0.17 


3 


0.7 


0.4 


DSK 


300 


CopB 


0000 


296 













UNC 


333 


CopB 


0000 


Note. 
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portion is shown here for guidance regarding its form and content. Column 1: X-ray source 
number. Column 2: SED slope from IRAC photometry with la error. Column 3: Number of 
IRAC bands from which the SED slope was derived. Columns 4-5: J-band source extinction and 
stellar mass estimated from the 2MASS J vs. J — H color-magnitude diagram and a 3(2) Myr PMS 
isochrones for the Cep OB3b/S 155(Cep B) regions assu r ning d — 725 pc. Column 6: Membership 
and PMS class. The membership is from fCetman et al 1 <2006l ') and the PMS class is derived here: 
"EXG" — possible extragalactic contaminant; "FRG" — possible foreground contaminant; "DSK" 
and "NoD" — disk-bearing and diskless PMS stars, respectively; "UNC" — object of the uncertain 
class. Column 7: Flag indicating transition disks as based on our visual inspection of SEDs. 
Column 8: Three-digit flag indicating source position on IR color-color diagrams: 3.6 — 4.5 vs. 
5.8 — 8.0 diagram, 3.6 — 4.5 vs. 4.5 — 5.8 diagram, and Ks — 3.6 vs. 3.6 — 4.5 diagram. Flag 
values: "0" — the PMS classification from Column 6 is consistent with location in the color- 
color diagram; "1" — disk-bearing PMS star from Column 6 lies in the locus of diskless stars in 
the color-color diagram; "2" — diskless PMS candidate from Column 6 lies in the locus of disk- 
bearing stars in the color-color diagram; "3" — not a PMS member from Column 6 or not enough 
information to place a member on the corresponding color-color diagram. Column 9: Subregion: 
"CepOB3b" — lightly absorbed Cepheus OB3b cluster; "CepB" — molecular cloud; "S155" — H II 
region interface. Column 10: Positional flag for each IRAC channel: "0" — source lies within the 
IRAC channel field of view; "1" — source lies outside the IRAC channel field. 



Table 3. Infrared Properties of non-Chandra IR-Excess Members 



No. 


R.A. 


Decl. 




J 


H 






Fi 




[3.6] 




[4.5] 


[5.8] 




[8.0] 


Subregion 


Aj 


M 




(deg) 


(deg) 




(mag) 


(mag) 




(mag) 






(mag) 




(mag) 


(mag) 




(mag) 




(mag) 


(Mq) 


(1) 


(2) 


(3) 




(4) 


(5) 




(6) 


(7) 




(8) 




(9) 


(10) 




(11) 


(12) 


(13) 


(14) 


1 


343.90431 


62.619685 


16, 


.26 ± 0.12 


14.99 ± 0. 


,09 


13.97 ± 0.06 


BAAccO 


12, 


,49 ± 0.02 


11 


.81 ± 0.02 


11.15 ± 0, 


,03 


10.34 ± 0.03 


CepOB3b 


1.3 


0.20 


2 


343.90796 


62.613150 


15, 


,67 ± 0.07 


14.36 ± 0, 


,05 


13.60 ± 0.04 


AAAOOO 


12, 


41 ± 0.02 


11 


.95 ± 0.02 


11.64 ± 0, 


,03 


11.06 ± 0.03 


CcpOB3b 


1.3 


0.30 


3 


343.90916 


62.601677 


15, 


,06 ± 0.04 


13.76 ± 0, 


,04 


12.85 ± 0.02 


AAAOOO 


11 


,66 ± 0.02 


11 


.12 ± 0.02 


10.67 ± 0, 


,03 


9.92 ± 0.03 


CcpOBSb 


1.2 


0.43 


4 


343.91215 


62.631087 


15, 


,20 ± 0.05 


14.21 ± 0, 


OS 


13.93 ± 0.05 


AAAOOO 


13, 


,32 ± 0.03 


12 


.98 ± 0.03 


12.57 ± 0, 


()3 


11.81 ± 0.03 


CcpOBSb 


0.7 


0.27 


5 


343.91421 


62.627495 


14, 


()8 ± 0.03 


12.97 ± 


,03 


12.29 ± 0.02 


AAAOOO 


11, 


,46 ± 0.02 


11 


.02 ± 0.02 


10.74 ± 


()3 


10.03 ± 0.03 


CcpOBSb 


0.7 


0.61 


6 


343.91472 


62.634968 


15, 


,22 ± 0.05 


14.20 ± 0, 


,05 


13.51 ± 0.04 


AAAOOO 


12, 


,66 ± 0.02 


12 


.29 ± 0.02 


12.14 ± 0, 


,03 


11.69 ± 0.03 


CepOB3b 


0.7 


0.27 


7 


343.92033 


62.641813 


15. 


.76 ± 0.07 


14.58 ± 0. 


.07 


13.88 ± 0.05 


AAAOOO 


12. 


.69 ± 0.02 


12, 


.28 ± 0.02 


11.95 ± 0. 


.03 


11.39 ± 0.03 


CepOBSb 


1.0 


0.24 


8 


343.92861 


62.671552 


15. 


.32 ± 0.05 


14.41 ± 0. 


.06 


14.06 ± 0.05 


AAAOOO 


13. 


.68 ± 0.02 


13, 


.36 ± 0.02 


13.22 ± 0. 


.03 


12.76 ± 0.04 


CepOBSb 


0.6 


0.22 


9 


343.93061 


62.645682 


16. 


.03 ± 0.10 


14.74 ± 0. 


.07 


13.88 ± 0.06 


AAAccO 


12. 


.67 ± 0.02 


12, 


.04 ± 0.02 


11.38 ± 0. 


.03 


10.36 ± 0.03 


CepOBSb 


l.S 


0.24 


10 


343.93414 


62.631564 


15. 


.21 ± 0.05 


14.24 ± 0. 


.06 


13.46 ± 0.03 


AAAOOO 


12. 


.18 ± 0.02 


11 


.72 ± 0.02 


11.48 ± 0. 


.03 


11.00 ± 0.03 


CepOBSb 


0.6 


0.25 



00 



Note. — This tabic is available in its entirety in the electronic edition of the journal. A portion is shown here for guidance regarding its form and content. Column 1: IR-cxccss source number. 
Columns 2-3: IRAC right ascension and declination for epoch ,12000.0 in degrees. Columns 4-7: 2MASS J H K ^. magnitudes, and 2MASS photometry quality and confusion-contamination flag. 
Columns 8-11: IRAC magnitudes derived in this work. Column 12: Indicates on which region the source position is projected: "CepOBSb" — lightly absorbed Cepheus OB3b cluster; "CepB" — 
molecular cloud; "S155" — H II interface. Columns 13-14: J-band source extinction and stellar mass derived from the J vs. J — H color- magnitude diagram by dereddening 2MASS photometric 
colors to the 3(2) Myr PMS isochrones for the Cep OB3b/S 155(CepB) regions assuming d = 725 pc. 



Table 4. Disk Fractions of mass complete PMS samples 



Class 


Ccp OB3b 


Ccp OB3b 


S 155 




S 155 


Cep B 


Cep B 




NE 


SW 


NE 




SW 


NE 


SW 






W) 






(5) 




\') 








0.5 < M < 1 M0 










DSK 


28 


14 


^ 10 




9 


5 


4 


NoD 


34 


14 


14 




4 


4 


3 


Disk Frac. 


0.46+°°^ 




42 + "-^" 


0. 


fiO+O.lO 
""-0.14 




o-57lS:'il 








M> 1 Mq 










DSK 


10 


11 


4 




3 


9 


4 


NoD 


16 


13 






3 






Disk Frac. 


n o„+0.10 


46+°-^ 


0.80+°0^ 


0. 


.,, + 0,18 
-0,17 




o.8o+°:5| 






Total over 


the M > 0.5 Mq 


mass range 






DSK 


38 


25 


14 




12 


14 


8 


NoD 


50 




15 




7 


6 


4 


Disk Frac. 


„ 43+0.05 


0.48+°;S^, 


n 4R+0 09 


0. 


RO+O.IO 

°-^-0.12 


'"-0.11 


"•°'-0.15 


Class 


Inner 




Intermediate 






Outer 






( 0'- 3') 




(4'- 7') 






(> 8') 










0.5 < M < 1 Mq 










DSK 


9 




~ 30 






10 




NoD 


7 




24 






20 




Disk Frac. 


„ _„+0.11 




















M> 1 Mq 










DSK 


13 




17 






2 




NoD 






13 






13 




Disk Frac. 


O.Sl+S.fs 




„ __+0.08 
"■°'-0.09 






„ ,,+0.13 
"■^■5-0. 04 








Total over 


the M > 0.5 Mq 


mass range 






DSK 


22 




47 






12 




NoD 


10 




37 






33 




Disk Frac. 


69+°'°'' 










"•^'-0.05 
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Fig. 1. — (a) Large-scale optical image from the Digitized Sky Survey (DSS) of the 
Cepheus B/0B 3 environs. Magenta contours and core labels outline ^^CO emission from 



SargentI (119771 ). O stars (cyan circles), B0-B3 (green), B4-B9 (red) are from iBlaauw et al. 



The black dashed line demarcates the older sub-group Cep 0B3a from Cep 0B3b. 
North is up, east is to the left, (b) Closeup of the Cep B region showing the 20' x 15' IRAC 
image in the 4.5 fim band. The ionizing sources of the region, the 07 star HD 217086 and 
Bl star HP 2170 6 1, are marked. The magenta contours show the ^^CO emission measured 
by iBeuther et al.l (120001 ) with the central hot core marked. The 17' x 17' Chandra ACIS-I 
field is outlined in cyan. The H II S 155 interface between the Cep B molecular cloud and 
unobscured Cep 0B3b association is schematically indicated in cyan. 
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Fig. 2. — Combined X-ray and IR image of the Cep B/Cep 0B3b region. The adaptively 
smoothed Chandra ACIS-I image in the 0.5 — 8.0 keV band (blue) is superposed on the 
Spitzer IRAC composite image in the 3.6 /xm (green) and 5.8 fim (red) bands. The 17' x 17' 
Chandra ACIS-I field is outlined in blue and the ~ 20' x 15' field of the Spitzer IRAC 
mosaic in 3.6/5.8 fim bands is outlined in green. 
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Fig. 3. — IR spectral energy distributions for nine representative diskfree (Class III) X-ray 
stars. From top to bottom, SEDs for stars in Cep 0B3b, S 155 and Cep B subregions. JHKg 
(diamond) and IRAC-band (square) flux points with usually small error s. The dash e d and 
solid lines give the original and (de)reddened IC 348 median SED from iLada et al.l (120061 ) 



fitted to the data. The top two lines of the panel labels give information from Table 



third line gives the spectral class of the IC 348 median SED template from lLada et al. 



a The 



(120061 ) 



and reddening applied to the original template SED to fit the observed Cepheus source SED. 
See the electronic edition of the Journal for the full atlas of IR SEDs of the Cepheus X-ray 
young stars. 
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Fig. 4. — IR spectral energy distributions for nine representative disk-bearing (Class II) 
X-ray stars. See Figure [3] for details. 
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Fig. 5. — IR spectral energy distributions for nine representative transition disk (Class II/III) 
X-ray stars. See Figure [3] for details. 
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Fig. 6. — Histograms of observed SED spectral slopes ao for (a) X-ray stars in the Cepheus 
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indicate number of stars with available photometry. 
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Fig. 8. — IRAC color-color diagrams of non- X-ray IR-excess sources. Top: [3.6] — [4.5] versus 
[5.8] — [8.0] color-color diagrams employed to select PMS candidates for (a) Cep 0B3b, (b) 
S 155, and (c) Cep B subregions. Dashed lines show classification criteria. Bottom: [4.5] 
versus [4.5] — [8.0] color-magnitude diagrams for (d) Cep 0B3b, (e) S 155, and (f) Cep B 
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vectors represent Ay ^ SO mag. 
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Fig. 9. — Mass distribution of Class III Chandra (red) and combined Class II Chandra and 
non-Chandra (blue) PMS stars in the (a) Cep 0B3b, (b) S 155, and (c) CepB sub-regions. 
The numbers of stars are estimated within mass bins of 0.2 dex width and la equivalent 
Poisson errors are shown. Annotation labels indicate the corresponding to each bin number 
of detected Chandra Class III stars (red), as well as Chandra + non-Chandra Class II stars 
(blue ). The dashed lines represent the scaled versions of the Galactic field IMF ( iKroupa 
20021 ). 
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Fig. 10. — Comparison of disk fraction as a function of stellar absolute J-band magnitude 
for four 2 — 3 Myr stellar clusters: Cep 0B3b (circles with la error bars); Cha 1 (solid 
line); IC 348 (dashed line); and a Ori (dashed-dotted line). The Cep 0B3b stellar sample is 
the combined Class 11/111 Chandra and C lass 11 non-Chandra P MS stars from this work. 
Data for the three other clusters are from iLuhman et al.l (j2008bl ). The top axis shows an 
approximate scaling to stellar mass. 
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Fig. 11. — Spatial distribution of Cepheus Chandra plus non-Chandra mass-complete se- 
lected PMS samples superposed on a grey-scale 5.8 /im band Spitzer-IRAC image: diskfree 
stars (red), disk-bearing stars (blue). Panel (a) shows 0.5 < M < 1.0 M0 stars, and panel 
(b) shows M > 1.0 Mq stars. The 07Vn star HD 2170786 is marked by the cyan square- 
circle, and the molecular hot core by the cyan cross. The cyan dashed locus delineates the 
S 155 subregion. The Chandra field of view is outhned by the large cyan square, and the 
NE and SW azimuthal zones are demarcated by the magenta line. 
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Fig. 12. — Running average of the disk fraction as a function of angular distance from 
Cep B hot core for the combined Chandra plus non-Chandra sample of PMS stars. The 
three subregions are separated by dashed lines and the green vertical arrow indicates the 
projected location of the primary ionizing star HD 217086. Mass-stratified sub-samples are 
color coded: M > 0.5 M© (red), 0.5 < M < 1 M© (black), and M > 1 Mq (blue). An 
interpretation based on three constant disk fraction layers is indicated by horizontal grey 
bars with their disk fraction levels consistent with those from Table HI 
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Fig. 13. — Vic color-magnitude diagram showing the Cep B/Cep 0B3b X-ray selected PMS 
stars individually corrected for extinction. Disk-bearing stars are plotted as open squares, 
and diskfree stars as filled circles. Colors indicate subregion location: Cep 0B3b (red), 
S 155 (green), and Cep B (blue). Solid lines show PMS isochro nes for ages 1, 2, 3 , 4, 5 Myr 
and the Zero- Age Main Sequence from evolutionary tracks of ISiess et al.l ( l2000l ) assuming 
d = 725 pc. Dashed lines show mass tracks 0.3, 1.0, and 3.0 Mq. The reddening vector 
shows Ay ~ 2 mag. 
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Fig. 14. — X-ray luminosity as a function of mas s for lightly-ob s cured PMS stars in (a) 
Cep 0B3b and (b) the Orion Nebula Cluster from iGetman et al.l (120051). In this figure, x 
symbols are diskless or non-accreting stars, • are optically thick disk or accreting stars, and 
□ are transition disk objects. Running medians for logL^, are shown as dashed (solid) lines 
for diskless or non-accreting (optically thick disk or accreting) stars. 



